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Abstract 



Mass splitting between axial and vector Qq mesons is considered within the standard 
QCD sum rules. In agreement with the first experimental data on the B\ meson (J p = 1 + ) 
we find that the splitting for B is about the same as for D and show that I/tuq corrections 
to the meson masses are small. 



Several groups [|TJ recently announced that they had observed candidates for the axial open 
beauty state (J? = lt, 2 ) about 500 MeV above the corresponding vector meson _B*(5325). 
The splitting appears to be of the same magnitude as in the case of D mesons, as predicted 
e.g. in the instanton liquid model||. This fact however looks somewhat strange from the 
viewpoint of expansion of the meson mass with respect to l/m,Q, where mq is the heavy 
quark mass, 

E l 

m = m Q + E -\ \- (1) 

m Q 

and indicates that the 1/tuq corrections are either very close in both the vector and axial 
cases, or small. On one hand, there appears to be no reason for these corrections to be 
close (e.g., the resonance energy, E , is different in both cases|J). On the other hand, it is 
known that for the couplings, such as fp and these corrections are very important, and 
1 / mg-expansion breaks down for /d @, • In the present paper we calculate the mass of the 
axial B\ meson and demonstrate, using a non-relativistic version 0, ||] of the standard QCD 
sum rules||, that I/tuq corrections to the Qq vector and axial meson masses are indeed 
rather small. 

The sum rules for vector and axial Qq mesons are obtained by considering correlators of 
vector (Jp = Q^^q) and axial (j^ = Qj^q) currents 

C, u {q) = * J ^e^'(0|T{j M (x), j+(0)|0> (2) 

at Euclidean momentum — q 2 > 1 GeV 2 . We choose the tensor structure proportional to 
g^u, since for q^q v the lowest state in the axial case is the pseudoscalar meson, and take into 
account operators with dimension d < 5 omitting the gluon condensate whose contribution 
is negligible. After the standard Borel transformation the sum rule takes the form|7], || 
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^exp(-m|/M 2 ) 

g± 

±m Q (qq)L^exp(-m 2 Q /M 2 ) =F -^m 3 Q m 2 (qq) exp(-m 2 Q / M 2 ) . (3) 

The couplings g± are defined according to (0|j /J |l ± ) = (m^/ g±)e fl , where e M is the meson 
polarization vector. The function F(x) decribes the a s -corrections||, 

F(x) = — + 2Li 2 (x)+lnxln(l-x) + --—ln(-—) 
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4 — x — x 2 5 — x — 2x 2 

— m(l — x) — -, r— 7 -xlnx — 7 -7 -. (4) 

v ; (l-x) 2 (2 + x) (l-x)(2 + x) V ; 

The other quantities in Eq.(g) are: m 2 , = {qa al3 G a(} q) / '(qq) = 0.8 GeV 2 , (gg) = -(0.24) 3 
GeV 3 , L = ln(M/A)/ln(///A), A = 0.15 GeV, p, = 0.5 GeV. The meson mass is obtained 
from Eq.(|3|) by taking the logarithmic derivative with respect to M 2 . The only difference 
between the vector and axial sum rule (apart from the value of the continuum threshold, sq) 
is the sign of the terms with (gg). Neglecting the continuum and the anomalous dimension 
factor, it is easy to obtain the following estimate 



m 2 + -m 2 _ ^ A7T 2 m Q (qq} t _ 5m 2 m 2 Q 



m 2 , + m 2 _ 



M 4 I VIM- 



(5) 



In the limit mq — > oo the Borel parameter scales as M 2 ~ 2mqfi, where /j, — 1/r is the non- 
relativistic Borel parameter and r is the typical euclidean time over which the correlators 
change significant lyH]. Thus, one expects for the splitting between axial and vector Qq 
states at mq — > oo 



m + -m^-— -U-^) • (6) 



The r.h.s. of the above equation has a flat maximum at \i ~ 0.4 GeV which corresponds to 
m + — m_ pa 0.4 GeV. The question is whether b and c quarks are heavy enough for B and 
D mesons to satisfy Eq.(^|), i.e. whether 1/mQ-corrections to masses are really small. Our 
experience with the couplings |j] indicates that these corrections may be important. 

In Figs. 1 and 2 we show the results of numerical analysis of the sum rules of Eq.(TJ) 
in the corresponding working windows in M 2 for m c = 1.35 GeV, m& = 4.7 GeV and the 
optimal values of continuum thresholds Sq = 6, 8, 35 and 40 GeV 2 for D*, Di, B* and Bi, 
respectively 0, R |TD[. We see that the mass splitting is indeed the same for D and B: 



m+ - m _ = (500 ± 50) MeV . (7) 

For the masses, the sum rules give the values rriB 1 = (5.9 ± 0.15) GeV, mg« = (5.4 ± 0.15) 
GeV, mflj = (2.5±0.15) GeV and mj). = (2.0±0.15) GeV. The errors correspond to allowed 
variations of continuum threshold Sq. A better accuracy for the splitting is due to partial 
cancelation of continuum contributions in this case. For the couplings we get gs 1 = 22 ± 3, 
qb* = 27 ± 3, go 1 = 10.5 ± 1.5 and g D * — 9 ± 1.5 (the last two couplings were obtained 
before 0, |§). It is worth mentioning that a s -corrections are rather important: without them 
the splittings would be about 30% bigger. 
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Now, we will explicitly demonstrate that the 1/mg-correction to the meson mass in 
Eq.(|l]) is small by expanding the sum rule of Eq.(^) around the limit rriq — > oo. Following 
the same procedure as in||, we introduce the non-relativistic continuum threshold E c and 
Borel parameter /x, s = (mq + E c ) 2 , M 2 = 2rriqfi, and obtain 

F 2 (Ejfi)a™ d ±TT 2 m 2 (qq)/U4^ 

tlQ = 6jl- 



' F±{E C / /i)a^ od ± n 2 (qq)L 4 / 9 /6fi 3 =F n 2 m 2 (qq) / 96 /i 5 



#1= - ^o'-T^exp^o/^/C^-^i^^c/M)^ 

2 Oo7T 

o 

4|Y Z7 1 TP \ n rad /oA^^dl 



exp[(£ " # c )//i]£ c 4 [(£o - - (H/^/2)ar] ■ (9) 



2C 7T 2 

where Co determines the asymptotic behavior of the residue in the limit rriq — > oo, 

J 5 = %fi + — + •■■) • do) 

The coefficient Co itself is determined from the following sum rule 

C = exp(£ /V) {^F 1 {E c /^)a\ ad ± (qq)L^ T rn%(qq)/16^ . (11) 

Upper and lower signs in Eqs. (§) and (|TT| ) correspond to axial and vector cases, respectively 
Once E± is known, 1/mQ-correction to the coupling can be obtained from the sum rule 



E 2 + 2E X 48 

+ 
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C x = - AE + ^ L- exp ( Eo /v)^F 2 (E c /vW 2 ad 



+ ^- 2 exp[(E -E c )/^]EX ad ■ (12) 

The functions F\ and F 2 in Eqs.(^), (0), ( [TTj) and (O) are the standard functions describing 
quark loop and continuum contributions, 



Fx(x) = 1 - (1 + x + x 2 /2)e~ x , 

F 2 (x) = 1 - (1 + x + x 2 /2 + x 3 /6)e~ x . (13) 

The factors a\ ad and a 2 ad contain radiative corrections [], 

1 As in ref.||, the arguments of the logs in a\ a 2 correspond to the maxima of the integrands in the 
dispersion integrals. In the numerical analysis we put vtiq = m&. Thus, we do not actually go to the limit 
m,Q — > oo in the radiative corrections. For a rigorous treatment of radiative corrections in this limit within 
the heavy quark effective theory, see e.g.[[Tl). 
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a\ ad = l + a s (^2.35 + - log(m Q /4/i)) , 

a r 2 ad = l + a s (2.35 + - log(m Q /6/x)) . (14) 

Note, that in the leading order in I/tuq the sum rules for J = 1 and J = meson masses |§ 
are the same. The hyperfine splitting is contained in the 1/mg-corrections. 

The asymptotic coefficients E and Co were calculated in ref.Q, and the splitting Eq — 
Eq = 800 ± 200 MeV was obtained. Our results for Eq and Eq are presented in Fig. 3 for 
the non-relativistic continuum thresholds E c = 1.8 GeV and E c = 1 GeV in the axial and 
vector cases respectively. We obtained a smaller value, Eq — Eq = 600 ± 100 MeV. We were 
not able to trace the source of this difference, since the axial case was discussed very briefly 
in ref. 0. For Cq we then get Cq ~ 0.71 and Cq ~ 0.16 (Fig. 4) which agrees with ref.[Q 
0. Using these values in Eq.(|9|) we obtain the results for 1/mq mass corrections presented 
in Fig.5. We see that the stability of the sum rule for E\ and the accuracy is rather poor 
in the axial case. However, it is clear that the 1/mQ-corrections turn out small both on the 
scale of rrib and m c 0. Finally, using the values Ef = and E± = —0.12 GeV 2 we get from 
Eq.(p|) Ct « -2.6 GeV and Cf w -5.5 GeV (Fig.6). 

Thus, we have shown that the splittings between axial and vector mesons with open charm 
and beauty calculated from QCD sum rules are rather close and agree with experiment. By 
a non-relativistic expansion of the sum rules we checked that 1/mQ-corrections to meson 
masses are very small. This is in contrast with similar corrections to the couplings which 
are very important (in the vector case 1/ttlq expansion does not work for «7b*!). 

I am grateful to A.B. Kaidalov for asking a question which triggered this calculation. 
I am indebted to H. Leutwyler for the warm hospitality at the University Bern where this 
work has been done. This work was supported in part by Schweizerischer Nationalfonds and 
by the INTAS grant 93-0283. 



2 This value of C is in agreement with the estimate obtained in an earlier paper 02] . 

3 1/ niQ corrections to vector and pseudoscalar meson masses were calculated within the sum rules approach 
in Refs.[^3) and [T^] . Our value for E\ in the vector case is of the same sign as in|jl|, but four times smaller 
in the absolute value. On the other hand, it is two times smaler in the absolute value and of the opposite 
sign than m[l4||. The disagreement is disturbing and should be resolved. But we do not discuss it here, since 
the corrections are small in any case. 
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Figure Captions: 

• Fig. 1: Mass splitting between D\ and D* from ([3]) in GeV. 

• Fig. 2: Same as Fig. 1 for B\ and B*. 

• Fig. 3: E in GeV. 

• Fig. 4: C in GeV 3 . 

• Fig. 5: E 1 in GeV 2 . 

• Fig. 6: C x in GeV. 
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Mass splitting for D 
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